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Abstract 

In this paper we continue our descriptions of stock markets in terms of some non 
abehan operators which are used to describe the portfoho of the various traders and 
other observable quantities. After a first prototype model with only two traders, we 
discuss a more reahstic model of market with an arbitrary number of traders. For 
both models we find approximated solutions for the time evolution of the portfolio 
of each trader. In particular, for the more realistic model, we use the stochastic 
limit approach and a fixed point like approximation. 
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I Introduction 



In a recent paper, [T], we have discussed how a quantum mechanical framework can be 
used in the analysis of stock markets. The conservation of the total number of shares and 
of the total amount of cash in any closed marked, i.e. in a market which does not interact 
with the environment, as well as the discrete nature of the number of the shares and of the 
monetary unit, suggested the use of some typical tools of QM^o, i.e. of quantum mechanics 
for systems with infinite degrees of freedom, in this different context. In particular, we 
have shown that a second quantized vision of the stock market produces in a natural way 
a set of differential equations describing the time evolution of the portfolio of each trader 
of the market. These results are on the same line as those given in [2] and [3], as well as 
in [1] and references therein. We should also mention that the use of tools coming from 
physics for economical problems, or more generally for dealing with complex systems, is 
a well established procedure, for which we refer to [5]. 

The paper is organized as follows: 

in the next section we briefly review the results in [1]. 

In Section III, contrarily with what has been made in [1], where the price of the share 
P{t) has no real dynamics since it is replaced by its mean value, we introduce a prototype 
model with only two traders were we show how to keep into account the time evolution 
of the price P{t). We prove that many integrals of motion exist. The equations of motion 
are solved using a perturbative expansion, well known in QM^. 

In Section IV we consider another model, which generalize the previous one in the 
sense that it consists of traders with arbitrary N > 2, and for which we consider two 
different approximations: the stochastic limit, which is useful to analyze the equilibrium 
of the model, and what we call a fixed point-lixe (FPL) approximation, which we use to 
deduce the approximated time evolution of the portfolio of any given trader of the stock 
market associated to the model. 

Section V contains our conclusions and plans for the future, while we discuss in the first 
Appendix some delicate mathematical points and in Appendix B a general introduction 
to the stochastic limit. 
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II The genesis of the model 

In this section we review some results and ideas first introduced in [1] which have produced 
an interesting toy model of a stock market based on the following assumptions: 

1. Our market consists of L traders exchanging a single kind of share; 

2. the total number of shares, A^, is fixed in time; 

3. a trader can only interact with a single other trader: i.e. the traders feel only a 
two-body interaction; 

4. the traders can only buy or sell one share in any single transaction; 

5. the price of the share changes with discrete steps, multiples of a given monetary 
unit; 

6. when the tendency of the market to sell a share, i.e. the market supply, increases 
then the price of the share decreases; 

7. for our convenience the supply is expressed in term of natural numbers; 

8. to simplify the notation, we take the monetary unit equal to 1. 

We refer to [T] for the analysis of these conditions, which however, in our opinion, look 
quite natural and self-explanating. The formal hamiltonian of the model is the following 
operator: 



I 





where P = p^p and the following commutation rules are assumed: 



[ai, al] = [q, 4] = 6inl, \p,p^ = [o, o^] = I, 



(2.2) 
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while all the other commutators are zero. We further assume that pa = 0. Here the 
operators aj, p^, cj and are respectively the number, the price, the cash and the supply 
operators, [1]. The states of the market are 

^{n};{k};0;Mi ■ ) =< V{n};{k};0;M ^ ■ ^{n};{k};0;M >, (2.3) 

where {n} = ni, 77-2, ... , til, {k} = ki, k2, ■ ■ ■ ,kL and 

(n^^^l . . . (alY^icl)''^ ■ ■ ■ (c\)''^(0^^^('niW 



{p^y 



•= ^ n,\...n,\k,\...kJ Um. ^^-^^ 

Here cpo is the vacuum of the model: ajifo = Cjipo = pipo = OLpo = 0, for j = 1,2, . . . , L. 
Again we refer to [1] or to any quantum mechanical textbook, see [6] for instance, for 
further details on second quantization. 

The interpretation of the hamiltonian is a key point in our approach and has been 
discussed in details in [1]: just as an example, the presence of the term p in Hj im- 
plies that when the supply increases then the price must decrease. Moreover, because of 
alaj{ciCj)^ , trader Tj increases of one unit the number of shares in his portfolio but, at 
the same time, his cash decreases because of cf , that is it must decrease of as many units 
of cash as the price operator P demands. Clearly, trader tj behaves in the opposite way: 
he loses one share because of aj but his cash increases because of (c])^. 

However we have discussed in [1] that the hamiltonian in (12. ip suffers of a technical 
problem: since Cj and cj are not self-adjoint operators, it was not obvious how to define 
the operators cj' and (c])^, and for this reason we have replaced H with an effective 
hamiltonian, H, defined as 

Hq + Hi, where 

= Zlf=i 0!ia]ai + Y,f=i I3ic\ci + o + p'^ p (2.5) 
= J2tj=iPij («l«i(ci c])^^ + ai a](cj clY'^^ + {o^ p + p^ o), 

where M = Lj^n};{k};0;M{P)- In this way, however, we are essentially freezing the price 
of our action, removing one of the (essential) degrees of freedom out from our market. 
This strong limitation will be removed in the next two sections of this paper and, in our 
opinion, this is really a major improvement. 
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Three integrals of motion for our model trivially exist: 




(2.6) 



This can be easily checked since the canonical commutation relations in fl2.2p imply that 



market. Of course, also the total amount of money must be a constant of motion since the 
cash is assumed to be used only to buy shares. Since also F commutes with moreover, 
if the mean value of o^o increases with time then necessarily the mean value of the price 
operator P = p^p must decrease and vice-versa. This is exactly the mechanism assumed in 
point 6. at the beginning of this section. Moreover, also the following operators commute 
with H and, as a consequence, are constant in time: 



The hamiltonian (12. 5p contains a contribution, hpo = o'^ o + p'^ p+ {o^ p + p^ o), which is 
decoupled from the other terms. For this reason it is easy to deduce the time dependence 
of both the price and the supply operators, as well as of their mean values. We get, [1], 



where we have called 0{t) = u;{n};{fc};0;A/(o^(i)o(^)) and Pr{t) = u;{n};{fc};0;A/(p'^(^)p(^))- 
Recall that Pr = Pr{0) = M. Equations ([23]) show that, if O = then 0{t) = Pr{t) = O 
for all t while, if O ^ then 0{t) and Pr(t) are almost constant. In the following we will 
replace Pr{t) with an integer value, the value M which appears in the hamiltonian (12.51) . 
which is therefore fixed after the solution (12. 8p is found. This value is obtained by taking 
a suitable mean of Pr(t) or working in one of the following assumptions: (i) O = Pr] or 
(ii) O ~ Pr or yet (iii) |0 + Pr| ^ \Pr — 0\. In these last two situations we may replace 
Pr{t), with a temporal mean, < Pr{t) >, since there is not much difference between these 
two quantities. However, in this way we are essentially removing the dynamics of the 




Qj — a] ttj ~^ 



(2.7) 



forj = l,2,...,L. 



P^{t) = 1{P,, + O + (P, - O) cos(2t)} 
0(t) = i{P, + O - {Pr - O) cos(2t)}. 



(2.8) 
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price from the model: no price variation occurs within this model after the replacement 
Pr(t) M\ How already anticipated, this restriction will be removed in the next section. 

Our main result in [1] was to deduct the time evolution of the portfolio operator, which 
we have defined as 

tij{t)=-inj{t) + kj{t). (2.9) 

Here we have introduced the value of the share 7 as decided by the market, which does 
not necessarily coincides with the amount of money which is payed to buy the share. As 
it is clear, Ilj{t) is the sum of the complete value of the shares, plus the cash. The fact 
that for each j the operator Qj is an integral of motion allows us to rewrite the operator 
tlj{t) only in terms of nj{t) and of the initial conditions. We find: 

Il,{t) = IljiO) + (7 - M){n^{t) - n,{0)), (2.10) 

In order to get the time behavior of the portfolio, therefore, it is enough to obtain fijit). 
We refer to [1] for a simple perturbative expansion of Hj(t) for L = 2. Here we prefer to 
show the other results, also contained in |T] , concerning the semiclassical thermodynamical 
limit of the model, i.e. a suitable limit for L 00. 

Our model is defined by the same hamiltonian as in (12. 5p but with M = 1. This is 
not a major requirement here since it corresponds to a renormalization of the price of 
the share, which we take equal to 1: if you buy a share, then your liquidity decreases of 
one unit while it increases, again of one unit, if you sell a share. Needless to say, this is 
strongly related to the fact that the original time-dependent price operator P{t) has been 
replaced by a certain weak mean value, M. 

It is clear that all the same integrals of motion as before exist: iV, i^, f , A := o — p 
and Qj = fij + kj, j = 1, 2, . . . , L. They all commute with H, which we now write as 

{H = h + hpo, where 
h = Yi=i Oiini + Y.Li Pi^i + E5=i («l«iCi c] + ai a]cj c| j (2.11) 
hpo = o"^ o + p^ p+ (o^ p + o). 

For hpo we can repeat the same argument as above and an explicit solution can be found 
which is completely independent of h. In particular we have 0J[n}\{ky,O;M{P) = 1- For this 
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reason, from now on, we will identify H only with h in 02.111) and we will work only with 
this hamiltonian. Let us introduce the operators 



X, = a,cl (2.12) 



1, 2, . . . , L. The hamiltonian h can be rewritten as 

L L 



1=1 i,j=l 

The following commutation relations hold: 

[X„ Xj] = 6,,{ki - n,), [X„ n,] = % X, [X,, k,] = -6,, X^, (2.14) 

which show how the operators {{Xj, Xj, fii, ki}, i = 1, 2, . . . , L} are closed under com- 
mutation relations. This is quite important, since, introducing the operators X^''^^ = 
'l2i=iPiiXi, / = 1, 2, . . . , L, we get the following system of differential equations, see [T]: 

Xi = t{f3i - ai)Xi + 2^x/'')(2ri, - Qi) 

4 = 2.(x,x/^)^-X(^)xt) ^'-''^ 

This system, as / takes all the values 1, 2, . . . , L, is a closed system of differential equations 
for which an unique solution surely exists. Indeed, we have found such a solution in [T] by 
introducing the so-called mean-field approximation which essentially consists in replacing 
Pij with ^, with p > 0. After this replacement we have that 



L 

1=1 i=l 

whose limit, for L diverging, only exists in suitable topologies, [TJ Ej, like, for instance, 
the strong one restricted to a set of relevant states. Let r be such a topology. We define 

_ L 

X°° = r - lim ^y^^i^ (2-16) 

i=l 

where, as it is clear, the dependence on the index I is lost because of the replacement 
Pii ^- This is a typical behavior of transactionally invariant quantum systems, where 
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Pi,i = Pi-i- The operator X°° commutes with all the elements of 21, the algebra of the 

observables of our stock market: = for all A G 21. In this limit system fl2.15p 

above becomes 

Xi = i{Pi-ai)Xi + 2iX°-{2n,-Qi) 

fii = 2i (Xi - X/) ^ ■ 
This system has been solved in under the hypothesis that 

Pi-ai=:<^^u (2.18) 

for all / = 1,2, ... ,L (but also in other and more general situations). We again refer to 
[1] for the details. Here we just write the final result, which is 

ni{t) = -1 - i^Y - 8\X^\^ iki{cos{ut) - 1) - ni{cosiut) + 1))} , (2.19) 

where we have introduced u = -\/($ — u)"^ + 16|X^|2. This allows also to find the time 
evolution for the portfolio, since Hiit) = 11/(0) + (7 — l){ni(t) — ri;(0)). Again, we refer to 
[1] for further comments and results. Here we just want to stress that our point of view 
has really produced, as an output, the time evolution of the portfolio of each trader of 
the market, which was indeed our original aim. 



Ill A two trader model 

As we have already discussed the model analyzed in the previous section has a very strong 
limitation: the time evolution of the price of the share, even if formally appears in the 
hamiltonian of the system, is frozen in order to get a well defined energy operator (i.e. 
in moving from H to H). Therefore, and in particular when we consider the thermody- 
namical limit of the model, such a dynamical behavior of the price operator completely 
disappears! 

In this section we will cure this anomaly, and for that we will discuss in many details 
a model based essentially on the same assumptions listed at the begin of Section II but in 
which the market consists of only two traders, Ti and T2. Of course, more than a realistic 
stock market, this can be seen as a sort of two-components physical system (ri + T2) 
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changing two different kind of particles (the shares and the money) and subjected to an 
external control (the price of the share and the supply of the system itself). However, 
even in view of the generalization which we will discuss in the next section, we will still 
refer to this physical system as a (toy model of a) stock market. 
The hamiltonian looks very much as the one in fl2.ip : 



H = Hq + Hi, where 



Hi = ^a\a2cf cl' + ai alc\ ' cl"^ + {o'^ p + p^ o) 



(3.1) 



with the standard commutation relations 

[o, o^] = \p,p^] = 1, [tti, a]] = [ai, a]] = 6ijl, (3.2) 

while all the other commutators are zero. 

The states of the system are defined as in (12.31) and (12.41) with L = 2, and the vectors 
<^{n};{fc};0;Af are cigeustates of the operators Ui = a|aj, ki = c\ci, i = 1,2, P = p^p and 
Cl = 0^0, respectively with eigenvalues ni,ki, i = 1,2, M and O. The main achievement 

here is that, how we discuss in Appendix A, we are now able to give a rigorous meaning 

,p 

to the operators Cj" and Cj , and this allow us not to replace the price operator with its 
mean value M and, as a consequence, to consider the price of the share as a real degree 
of freedom of the model. However, before defining cJ" and Cj , it is worth noticing that 
the non abelianity of our structure does not automatically implies that the observables 
of the market, i.e. the operators ki, hi, P and tl, as well as some of their combinations, 
cannot be measured simultaneously. This is because these observables, which are the only 
relevant variables for us, do commute and, as a consequence, they admit a common set 
of eigenstates, see equation (13. 3p below. 

Using the same arguments given in Appendix A we are able to define the operators 
Cj' and Cj via their action on the orthonormal (o.n.) basis of the Fock-Hilbert space Ti. 
of the model whose generic vector is, in analogy with (12.41) . 



\py 



^n.,nrMM;0;M - n,\ k,\ k,\ 0\ M ^^'^^ 



9 



Here rij, kj, O and M are non negative integers, (fo is the vacuum of the model, ajifo = 
Cjipo = pipo = oipo = 0, for j = 1,2, and H is the closure of the linear span of all 
these vectors. Then we have, for instance, ai</?ni,n2;fci,fc2;0;M = ^/nl'^n^-l,n2■,k^M^o■,M if 

ril > and ai(pnx,n2\kx,k2\0;M = if 77-1 =0, a\(pni^n2\kx,k2\0; M = \/ Ui + 1 (Pni+l,n2; ki,k2; O; M ■ 

Analogous expressions for the action of 02, al, Cj, c^j, o, o\ p and on (Pni,n2;ki,k2;0;M 
can also be recovered, see [6]. Moreover we have, see Appendix A, 

V5ni,n2;fci,fc2;0;M, if M = 0, V/^i > 0; 

c[^n,,n2;k,,k2;0;M:= { 0, if M > k,, ^k, > 0; (3.4) 



ki <^ni,n2;ki-M,k2;0;M, if /Cl > M > 



and 



h Vn^,n2;k-i_M\0;M '■— \ /TTTm y'^-^l 



where we have defined 



</2ni,n2;fci,fc2;0;M, if M = 0, Vfci > 0; 

k\^^^^ LPnt,n2;ki+MM\0\M-i if M > 0, 

k\~^'^ := ki{ki - 1) ■ ■ ■ (fci - M + 1) 
k\+''^ := (A;i + l)(A;i + 2) ■ ■ ■ {ki + M) 



(3.6) 



Analogous formulas hold for and cl . These definitions have a clear economical inter- 
pretation: acting with cf on (fnun2;kuk2;0;M returns (fnj_,n2;kuk2;0;M itself when M = 
since, in this case, the action of of coincides with that of the identity operator: the price 
of the share is zero so you don't need to pay for it and hence your cash does not change! 
Moreover, ii M > ki, destroys more quanta of money than ti really possesses. There- 
fore, the result of its action on the vector is zero. A similar problem does not occur when 

t ^ 

we consider the action of c\ on (Pni,n2;ki,k2;0]M, since in this case the cash is created! 

In the rest of the section, however, these formulas will be significantly simplified by 
assuming that, as it is reasonable, during the transactions between ri and T2 the price of 
the share never reach the zero value and, moreover, that no trader try to buy a share if 
he has not enough money to pay for it. Therefore we simply rewrite fl3.4p and 03.51) as 



Cl V'ni,n2;fci,A:2;0;M — Y fi^l '^ni,n2; ki-MM] O; M , fo 7\ 

'^ni,n2;ki,k2;0]M — V "'1 V^ni,n2; fci+M,fc2; O; M 
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The commutation rules are the standard ones, see fl2.2l) . plus the ones which extend the 
rules in (1A.4[) : 

[Ac/] = [P,c]V0, (3.8) 

and 

[k^cf] = -6,,Pcf, [fc^cf] = 6,,iPcf (3.9) 

for j = 1,2. 

Since our market is closed it is not surprising that the total number of shares and 
the total amount of cash are preserved. This is indeed proved simply computing the 
commutators of the total number of shares and the total cash operators, N = ni + n2 and 
K = ki + k2, with the hamiltonian H. Indeed one can check that [H, K] = [H, N] = 0. 
Moreover, we can also check that T := Cl + P commutes with the hamiltonian. This is, 
as already discussed in the previous section, the mechanism which fixes the price of the 
share within our simplified market: the more the market supply increases the less is the 
value of the share, i.e. its price. 

As already stressed before, one big difference between the model we are considering 
here and the one considered in Section II and in [1] is that now the price operator is not 
replaced by its mean value. This has an important consequence: the operators extending 
Qj in (12. 7p for this model, which are proportional to Pa'^aj + cjc-,-, j = 1,2, are no 
longer constants of motion, and they cannot be used to facilitate the computation of 
the portfolios of ti and T2. Nevertheless we will still be able to deduce, with an easy 
perturbative approach, the time behavior of both portfolios at least for small values of t. 

The first step consists in deducing the time evolution of the price of the share. This 
computation is completely analogous to that of jTj and will not be repeated here. Again, 
we can deduce that A := o — p is another constant of motion and we find that, see (12.81) . 

Pit) = uJn„nr,k^M;0;MipKt)pit)) = 1{M + Q + (M - O) cos(2t)} 

0{t) = Un„nr,kukr,O;M{0^{t)0{t)) = ^{M + O - [M - 0)cOs{2t)} 

In [1] the absence of a true dynamics for P suggested to define the portfolio of the j — th 
trader by introducing another parameter, 7, which was interpreted as the price of the 
share as decided by the market, which does not necessarily coincides with M. However, 
there was no direct link between 7 and M in [T], and this is not completely satisfying, of 
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course! Here we have no need for introducing such an extra parameter since we are now 
in a position to consider directly P{t) instead of its mean value M. Therefore we replace 
formula (12.91) by defining the portfolio of the trader tj as 

n,{t)=P{t)n,it) + k,it) (3.11) 

with j = 1, 2, which is just the sum of the total price of the shares and the cash of tj. 

Of course, due to the fact that P{t) is known, Ili{t) is known when we both know 
ni{t) and ki{t). Moreover, if we now ni{t) and ki{t), then we also know n2{t) and k2{t) 
since their sum must be constant, so that we can also find the analytic form of 112 (t). 
However, this is not the only way to find ni(t). Another possibility follows from the fact 
that, as it is easy to check, 

fli(t) = >(t)ni(t), (3.12) 

which shows again, even without any need of using Qj as in the previous section, that it 
is enough to know ni(t) to find the time evolution of the portfolio of Ti. 

However, even for this two-traders model, it is not easy to deduce the exact expression 
for Hi(t). Nevertheless, a lot of information can be deduced, mainly for short time 
behavior, using different perturbative strategies. Here we just consider the most direct 
technique, i.e. the following perturbative expansion 

Mt) = e*^*ni(0)e-^^* = fli(O) + zt[H, 11,(0)] + ^-^[H, [H, fli(O)]] + ■ ■ ■ , (3.13) 

leaving to the next section a more detailed analysis of other strategies to produce tlj{t). 
The computation of the various terms of this expansion, and of their mean values on 
the state uJni,n2;ki,k2;0;M{-), IS bascd on the commutation rules we have seen before and 
produce, up to the second order in t, the following result: 

ni(t) = UJ^,,nr,k,M;0;M{tll{t)) = Hi(0) + - M), 

which shows that, for sufficiently small values of t, the value of Ili{t) increases with time 
if O > M, i.e. if at t = and in our units the supply of the market is larger than the price 
of the share. It is further possible to check that the next term in the expansion above is 
proportional to An^^n2;ki,k2\M where 
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We avoid the details of this computation here since they are not very interesting, mainly 
because this is just a toy model which is more important for its general structure than for 
a real financial interpretation. Here we just want to stress that the expansion in fl3.13p 
gives, in principle, the expression of ni(t) at any desired approximation. 

IV Many traders 

In the previous section we have learned how to define the operators and its adjoint and 
we have used this definition in the analysis of a simple hamiltonian which was essentially 
already introduced in [1]. We devote this section to a more realistic model, where the 
stock market is made of different traders with arbitrarily large. 

In our approach we will focus our attention on a single trader, r, which interact with 
an ensemble of other traders in a way that extends the interaction introduced in (13.11) . In 
other words we divide the stock market, which as before is defined in terms of the number 
of a single type of shares, the cash, the price of the shares and the supply of the market, 
in two main ingredients: we call system, S, all the dynamical quantities which refer to a 
fixed trader r: its shares number operators, a, and n = a, the cash operators of r, 
c, and = c as well as the price operators of the shares, p, p^ and P = p^ p. On 
the other hand, we associate to the reservoir, TZ, all the other quantities, that is first of 
all, the shares number operators, A^, A|, and Nk = ^fc and the cash operators, Ck, Cl 
and Kk = Cl Ck of the other traders. Here k E A and A is a subset of N which labels 
the traders of the market (other than r). It is clear that the cardinality of A is A^ — 1. 
Moreover we associate to the reservoir also the supply of the market, which is described 
by the following operators o^, o|, and Ok = o^Ok, k E A. The stock market is given by 
the union of S and TZ, and the hamiltonian, which extends the one in (13. ip . is assumed 
here to be 

{H = Ho + \Hj, where 
Ho = Uan + u,k + ujpP + EfcgA {^Aik) Nk + ricik) Kk + ^o{k) Ofc) (4.1) 
Hj = (zt Z{f) + zZ\J)) + {p^o{g) +po^{g)) 

Here Ua, ojc and ojp are positive real numbers and VLA{k), VLc{k) and fio(^) are real valued 
non negative functions, whose interpretation was first discussed in |T]: they describe the 
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free time evolution of the different operators of the market. We have also introduced the 
following smeared fields of the reservoir: 

' ^(/) = EkeA Zk f{k)_= E.GA A, ci fiki 

< ZKf) = E,eAZlf{k) = E,eAAiCk^fik) (4.2) 

as well as the operators z = ac' , = CI. and their conjugates, since for instance 
Ak and Ck appear always in this combination both in Hj and in all the computations we 
will perform in the following. This is natural because of the physical meaning of, e.g., z: 
the action of 2; on a fixed vector number destroys a share in the portfolio of r and, at 
the same time, creates as many monetary units as P prescribes! Of course, in Hj such 
an operator is associated to Z'^^f) which acts exactly in the opposite way on the traders 
of the reservoir: one share is created in the cumulative portfolio of TZ while P quanta of 
money are destroyed, since they are used to pay for the share. The following non trivial 
commutation rules are assumed: 

[c, ct] = [p,pt] = [a, at] = 1^ [o,^^ ^t] = [A,, A]] = [Q, Cj] = 6,,1 (4.3) 
which implies 

[k,, C[] = -P Cf 4,., [kk, Cl ^] = P Cl ^ 4,, (4.4) 

Finally, the functions f{k) and g{k) in fl4.ip and fl4.2p are sufficiently regular to allows 
for the sums in 04.21) to be well defined, as well as the quantities which will be defined 
below, see (14.101) . 

Remark:— Of course, since r can be chosen arbitrarily, the asymmetry of the model is 
just apparent. In fact, changing r, we will be able, in principle, to find the time evolution 
of the portfolio of each trader of the stock market. 

The interpretation suggested above concerning z and Z{f) are also based on the 
following results: let 

N:=n + Y,Nk, K := k + t := P + Y.^k (4-5) 

fcGA fceA fcGA 
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Of course N is associated to the total number of shares in our closed market and, therefore, 
is called the total number operator. K is the total amount of money present in the market 
and is called the total cash operator. T has not a direct interpretation so far, since is 
just the sum of the price and the total supply operators, O = XlfceA^*^- worth 
recalling that the supply operators are only related to the reservoir TZ, because of our 
initial choice. This is the reason why there is no contribution to the operator O coming 
from S. 

Proposition 1 The operators N, K and f are constants of motion. 

The proof of this proposition is a simple exercise based on the commutation rules above. 
Indeed, it is not hard to check that H commutes with A^, K and with F. This proves 
that our main motivation for introducing the hamiltonian in (14.11) is correct: with this 
choice we are constructing a closed market in which the total amount of money and the 
total number of shares are preserved and in which, if the total supply increases, then the 
price of the share must decrease in order for F to stay constant. Of course, it would be 
interesting to relate the changes of O to other (maybe external) conditions, but this will 
problem will be considered elsewhere: here we just consider the simplified point of view 
for which O may change in time, but we don't analyze the reason why this happens. 

The next step of our analysis should be to recover the equations of motion for the 
portfolio of the trader r, defined in analogy with (13.111) as 

flit) = P{t)h{t) + kit). (4.6) 

It is not surprising that this cannot be done exactly so that some perturbative technique 
is needed. We will consider in the following sub-sections two orthogonal approaches, 
orthogonal in the sense that they give different information under different conditions 
which, together, help in a better understanding of the model. In particular we will first 
consider the so-called stochastic limit of the system: this approximation will produce the 
explicit form of the generator of the semigroup arising from the hamiltonian (14. ip . and 
this will give some interesting condition for the stationarity of the model, i.e. for n(t) 
to be constant in time. We will see that this is possible under certain conditions on the 
parameters defining the model. The second approach will make use of a sort of FPL 
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approximation which will produce a system of differential equation for the mean value of 
n(t) whose solution can be explicitly found. 

IV. 1 the stochastic hmit of the model 

The stochastic limit of a quantum system is a perturbative strategy widely discussed in 
[S] and which proved to be quite useful in the analysis of several quantum mechanical 
systems, see [10] for a recent review of some applications of this procedure to many-body 
systems. 

Here we adopt this procedure pragmatically, i.e. without discussing any detail, while, 
to keep the paper self-contained, we postpone to Appendix B the list of some basic facts 
of this approach. 

The first step consists in obtaining the free time evolution of the interaction hamilto- 
nian which we still call, with a small abuse of language, Hi{t). Due to the commutation 
rules (Oi) and ( l44l) we find that 

Hi{t) := e'^^'Hie-'"' = Z{fe''^^) + z Z^J e-''^^) + o{g e'''") + p o\g e"'''"), (4.7) 

where we have defined 

iz{k):=P{nc{k)-u,)-{nA{k)-LUa), eo{k) := cUp - no{k) (4.8) 

and, for instance, Z{f e'^^^) = EfceA/(^) ^'^'^''''^ ^k- 

The next step consists in computing first uj (if/ (p-) Hi (p-)), then 

and finally the limit of I\{t) for A ^ 0. Here a; is a state of the market, which we 
take as a product state oj = Usys ® i^res with Usys a gaussian state, that is it satisfies 
Usysia^) = uJsys{c^) = i^sysip^) = aud Usysiaa) = Usys{,cc) = Usysia^ a^) = Usysipp) = 
WsysiP^ P^) — 0- Here a" can be a or and the same notation is adopted for c" and 
pK These conditions are obviously satisfied if Usys is a vector state analogous to that in 
( 12.31) . We don't give here the details of the computation, which is rather straightforward, 
but only the final result which is obtained under the assumptions that the two functions 
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ez{k) := uj{ez{k)) and eo{k) are not identically zero. Moreover, it is convenient to assume 
that 

eY{k)=eY{q) ^ k = q, (4.9) 
where Y = Z,0. Then, if we define the following complex constants 



o 



(4.10) 



which surely exist if f{k) and 5f(/c) are regular enough, we get 

J(t) = -t ^^Usysiz^ z)T^z'^ + Usysiz Z^)T^z'^ + Usysip^ p)T''^^ + ^ sys{p )'^o^ } 



Next we need to find the expression of a self-adjoint, time dependent operator H^''^\t), 
the so-called stochastic limit hamiltonian, which reproduces this result in a sense that we 
will specify in a moment. 
Let us take 



H^^'\t) = (Z^''\t) + Z(^)^(t) ) + z (Z^''^\t) + Z^^\t)] + 



+p^ (o('^)(t) + o^'^V)) +p{o^''^\t) + o(')(t) 



(4.11) 



where the new operators introduced here are assumed to satisfy the following commutation 
rules: 

Z('^)(t), Z('^)^(t')l = Tf 5+{t - t'), \z^'\t), Z('')^(t')l = r^') 6+{t - t'), (4.12) 



and 



o('^)(^),o^"^'(0 =r'^'s+{t-t'), o^')(t),oW'(t') =T);;>64t-t') 



(4.13) 



if t > t'. The time ordering is crucial here and 5+ is essentially the Dirac delta functions 
but for a normalization which arises because of the time ordering we consider here, [9]. 
The only property of S+ which we will need is the following: 5+(t — r) /i(r) dr = h{t). 
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Now, let ^0 be the vacuum of the operators Z^^'^t), Z^^\t), o'^'^^t) and o^^\t). This 
means that Z('^)(t)^o = ^^^Hi)^o = o(")(t)*o = o^^\t)'^o = for all t > 0. Then, if we 
consider = Usysi-)'^ < ^o, • > and we compute 

Jo Jo 

we conclude that J{t) = I{t). This means that, at a first order, H^^^\t) allows us to get 
the same wave operator Ut which describes the time evolution of the systems. We use 
H^'''^\t) to construct the wave operator as f/^ = 1 — ? H^'''^\t') U^, and then to deduce 
the following commutation rules: 

[Z^''\t),Ut] =-zrfzUt, [Z^'\t),Ut] =-iTfz^Ut (4.14) 

and 

[o^''\t),Ut] = -zrg^f/t, [o^'\t),Ut] = -iT^Sv'Ut (4.15) 

by making use of the time consecutive principle, [9]. 

We are now ready to get the expression of the generator. Let X be a generic observable 
of the system, that is, in our present context, some dynamical variable related to the trader 
T. Let Ir be the identity operator of the reservoir. Then the time evolution oi X in 
the interaction picture is given by jt{X ® 1^,) = U}{X ® ^r)Ut, so that 

dtUX ® 1,) = zf//[if(^^)(t),X ® l,.]Ut 

Using now the commutators in fl4.14p and (I4.15p . and recalling that \l/o is annihilated by 
all the neui reservoir operators, we find that 



O {dtjtiX ® 1,)) = n{uHT^z\^\X] z - T^z^ [z,X] + Tf[z,X] z^ - vfz [z\X] + 



+T^S^ [pi x]p- T^S^p^ [p, X] + rg) [p, X] - rgV [pt, xm) 

which, together with the equality f2 ((9jjj(X (g) 1^)) = Q{jt{L{X ® Ir))), gives us the 
following expression of the generator: 



L{X ® 1,) = ^^"^[;z^X] z - T^z^z^ [z,X] + Tf[z,X] z^ - vfz [z\X] + 
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+ rg) [pt, x]p- rgV^ [P, X] + rg^ [p, x] pt _ rgV [pt, x]} (4.i6) 

Therefore we find, after few computations, 

L{n ® Ir) = 2^{Tf}zz^ -2^{T^^^}z^ z, (4.17) 

and 

L(^® 1,) = -2p^{r'^z^}zz^ + 2p^{r'^z^}z^ z, (4.18) 

which in particular shows that L{k (g) 1^) + P L{n ® 1^.) = 0. Finally we find, using these 
results and recalling that ti{t) = P{t) h{t) + k{t), 

L(fl® 1,) = 2(3?{rg^}-3?{rg^})Pn + 2 3?{rg^}n. (4.19) 

The first remark is that, in the stochastic limit, even if the time dependence of h and 
k depend on F^'' and F^'', the time evolution of fl in a first approximation do not! In 
fact, formula fl4.19p shows that it only depends on T^q and Fq\ However, since the 
time evolution of n depends on r|^'^^ because of (14.171) and (14.181) . this dependence will 
necessarily play a role also in n(t). 

The above equations show that, even after the stochastic limit has been taken, it is 
quite difficult to produce a closed set of differential equations. On the contrary it is quite 
easy to deduce conditions for the stationarity of the market. This is exactly what we will 
discuss next. 

We begin noticing that, for instance, we have 

2^{T^;}} = Y,\9{k)\^^res{oko{) I rfre-^^^o(^) = 27r5^|(7(A;)|2u;,es(ofeoi)5(£o(A;)) 

(4.20) 

and analogously we find that 3f?{Fg^} = T^Ylk&K \g{k)\'^ uJres{o\ok) 5{eo{k)) while ^{Vf} = 
7TEkeA\fik)\'ujUZ,Zl)6{ez{k)) and ^{vf} = TrZkeAlfi^)? ooUZl Z,) 6{ez{k)). 

Therefore, since [ok, ol] = 1 and, as a consequence, Uresiok ol)—uJres{ol Ok) = UresW = 
1, we find that ^{Tq^} — ^{Tq^} = —'n'J^keA l5'(^)P^(^o(^))- The conclusion now follows 
from (14.191) : the portfolio of r is stationary (in our approximation) when the function 
eo{k) has no zero for k E A. Indeed, if this is the case, we deduce that L(II ® 1^) = 0. 
Since eo{k) = ujp — Qo{k) this means that if the free dynamics of the price and the supply 
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are based on substantially different quantities then the portfolio of r keeps its original 
value, even if the operators n(t) and k(t) may separately change with time. This is an 
interesting result since it can be summarized just stating that, within the approximation 
we are considering here, the fact that n(t) depends or not on time is only related to a 
given equilibrium, if any, between the free price hamiltonian, ujpp^ and the free supply 
hamiltonian, XlfceA ^o(^) c^: again, the interplay between these two ingredients of the 
model play an interesting role! 

A similar analysis can be carried out also to get conditions for the equilibrium of h{t) 
and k{t). Because of (14.171) and fl4.18p . and because of the known time evolution of -P(t), 
h{t) is constant if and only if k{t) is constant, and for this to be true the function ez{k) 
must be different from zero for each /c G A. On the other hand, if at least one zero of 
ez{k) exists in A, then a non- stationary condition for n{t) and k{t) is possible. 

IV. 2 a different approximation 

The approach we have discussed so far produced some interesting information about the 
stationarity of the portfolio of r but no concrete insight about its time evolution. In other 
words, if we try to deduce the time behavior of 11 (t) we get no significant simplification 
if we adopt the form of the generator in (14.161) or if we look directly to the Heisemberg 
expression for n(t), n(t) = e*^* 11(0) e~*^*. However, also this last attempt does not 
produce directly a closed system of differential equations: some different approximation 
must be assumed. This different approximation will be discussed in this subsection. 

We first remind that given a generic operator X its time evolution, in the Heisemberg 
representation, is (formally) given by X{t) = e*^* X e"*''^*, and it satisfies the following 
Heisemberg equation of motion: X(t) = ie^^^[H, X]e~'^^^ = i[H,X(t)]. In the attempt 
of deducing the analytic expression for H(t), the following differential equations can be 
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deduced: 



^^ = ^X{-z^{t)Z{f,t) + z{t)Z^(f,t)), 
4f) = zXP{t) {z^{t) ZU\t) - z{t) Z\J,t)) , 
^^ = zA {p{t)oKg,t)-p^{t)o{g,t)), 

^ = I (P{t)uJ, - Ua) Z{t) + l\[z^{t), Z{t)] Z{f, t), 

dZ{U) 
dt 



I Z {Pmc - ^A)f,t) + a zit) [Zt(/, t), Z{f, t)]. 



where we have defined Z{f, t) := e'^' Z{f)e-'^\ Z ({P{t)Qc - ^A)f, tj = EfceA(^W^c(^) 
^Aik)) f{k) Zkit), o{g,t) = e'^^o{g) e"*^*, and so on. 

It is clear that the system fl4.2ip is not closed, since for instance the differential equa- 
tion for P{t) involves p{t), o{g, t) and their adjoint. This is not a major problem since, as 
in Sections II and III and in [1], it is quite easy to deduce the time evolution of the price 
operator P with no approximation at all. This is because p{t) (and o{g,t)) can be found 
explicitly. Even if these operators can be found under more general conditions, we will 
now restrict the model requiring that the coefficients in H satisfy some extra requirement, 
which are only useful to simplify the computations. For instance, we will assume that 
Qo{k) is constant in k, Qq = ^o{k) for all k E A, and that Up = J2keA Ifl'l^)!^ = A = Qq- 
Then we get p{t) = \ {p{e-^'^^ + 1) + o{g){e-^'^^ - 1)) and P{t) = p^{t)p{t). Since P{t) 
depends only on the operators p and o, and not on a, c, and so on, and since we are 
interested to the mean value of the operators in (14.211) in a vector state u generalizing 
(12.31) . we replace this system with its semiclassical approximation 



dh{t) 
dt 

dk{t) _ 

dt ~ 
dz{t) _ 

dt ~ 
dZ{f,t) 
dt 



t\{-z^it)Z{f,t) + z{t) Zt(7,t)), 
z\P,{t) {z^{t)Z{f,t)-z{t)Z^(f,t)), 

I (Pe(t)^c - UJa) Z{t) + l\[z^{t), z{t)] Z{f,t), 

= t Z - ^A)f. t) + z{t) [Zt(7, t), Z(/, t)]. 



(4.22) 



where 



P^{t) = uj{P{t)) = I [(M + O) + (M - O) cos(2At)] 



(4.23) 



We refer to P for a more complete discussion of the two-fold role of the state u. Here we 
just want to remark that a given vector state allows us to pass from the quantum dynamics 
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of the model to its classical counterpart, since we use uj to replace the time dependent 
operators with their mean values, which are functions of time. At the same time, moreover, 
a vector state is used to fix the initial conditions of the differential equations, that is the 
initial number of shares, the initial cash and so on. 

In order to simplify further the analysis of this system, it is also convenient to assume 
that both Qc{k) and fi^l^) are constant for A; G A. Indeed, under this assumption, the 
last two equation in fl4.22p forms by themselves a closed system of differential equations 
in the non abelian variables z{t) and Z{f,t): 



= I {P,mc - ^a) Zif, t) + z{t) [Zt(7, t), t)]. 



(4.24) 



Getting the exact solution of the system (14.221) . with (I4.24p as the two last equations, is 
an hard job. However, this is a good starting point for finding an approximated solution 
of the dynamical problem. Indeed, a natural approach consists in taking the first non 
trivial contribution of the system, as usually done in perturbation theory. This means 
that, in system (14.241) . the contributions containing the commutators must be neglected 
since they are proportional to A while i {Pc(t)uJc — uja) z{t) and i {Pc{t)Qc ~ ^a) Z{f,t) 
which, on the other way do not depend on A, give a relevant contribution. On the other 
way, in order not to trivialize the system, we have to keep the first two equations in (I4.22p 
as they are: if we simply put A = here, in fact, we would trivialize the time evolution 
of both n{t) and k{t). With this choice we get 

^ = ^X{-z1(t)Z{f,t) + zit)Z^{f,t))_, 
^ = .A Pe(t) (.t(t) Z{f, t) - z{t) Z^(f, t)) , 
^ = ^(P,(t)a;e-^a) zit), 

^ = 2 (Pe(t)l^C7 - f^A) Zif,t). 

However, we will now show that this approximation is too rude, meaning with this that, 
even if the operators n{t) and k{t) have a non trivial dynamics, at the classical level we 
deduce that both n{t) = uj{n{t)) and k{t) = uj{k{t)) are constant in time, so that the time 
behavior of the portfolio H(t) = Pc(t) n{t) + k{t) = Pdf) n + k is uniquely given by Pdt). 
We first observe that z{t) and Z{f,t) in (14.251) are 

zit) =ze'^^'\ Zif,t) = Z(/)e^*W, (4.26) 



22 



where 

X{t) = at + /3sin(2At), = at + Psm{2Xt) (4.27) 

with 

( a = |((M + 0)uj, - 2^J, P=t,{M- O) ^^ ^g) 



a = 1((M + 0)nc - ^Qa), P = ^(M - O) 

Our claim is now an immediate consequence of fl4.26p above. Indeed, from the first 
equation in (I4.25p . taking its mean value on the number vector state u we find 

h{t) = j^^m) = ^ (I^W) = {-^ {^Ht)Z{f,t)) + CO {z{t)Z^{f,t))} = 

since, for instance, u (^z^^ {t)Z{f,t)) = e~^^^^*^~^^*^^uj{z^ Z{f)) = 0. Analogously we find 
that k(t) = ^ uj(k{t)) = 0. Therefore we have n{t) = n and k(t) = k, as claimed above. 

Remark: in a certain sense this result relates the two approximations considered so 
far. Indeed, replacing (14.221) with (I4.25P we obtain a stationary behavior for n(t) and 
k{t). An analogous behavior was deduced, in the previous subsection, if ez{k) has no 
zero. However, these two different approximations cannot be directly compared. The 
reason is the following: in the stochastic limit approach we need to require that ez{k) 
and eo{k) are not identically zero. This is crucial to ensure the existence of lim^^o -^a(^)- 
In the present approximation we are requiring that both Qcik) and QA{k) are constant 
in k so that, see (14. Sp . we would get ez{k) = Pc{t){Qc — ^c) — (^a — cua), which may have 
some zero in k only if it is identically zero in k. In other words, in the conditions in which 
we are working here the stochastic limit approach does not work. Viceversa, if we are in 
the assumptions of the previous subsection, then system (14.220 cannot be easily solved! 
Hence the two approximations cover different situations. 

A better approximation can be constructed. Again the starting point is the system 
(I4.24p . for which we now construct iteratively a solution, stopping at the first relevant 
order. In other words, we take Zo{t) and Zq^/, t) as in (14.261) . Zo(t) = z e*^^*) and Z^^j , t) = 
Z(f)e^^^^\ and then we look for the next approximation by considering the following 
system: 

^ = Z {P,{t)uJ, - UJa) Zo{t) + l\ Zo(/, t) [zl{t)Mt)l 

= t{P,mc-nA)Zo{f) + t\zo{t) [Zl{f,t),Zo{f,t)]. 
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which can be still written as 

= I {P,{t)u, - Ua) Zo{t) + l\ Zoif, t) [zl Zol 



^( - V- cv/~c -aj ~KJ\-J i - ■ "UW ) L~U' -ujl f4 29) 

§^ = t{p,{t)nc-nA)Zoif) + tXzo{t)[zl(f),Zo{f)]. 



These equations can be solved and the solution can be written as 

z,it) = z viit) + Z{f) [z\ z] r/2(t), t) = Z{f) Ut) + z [Z(7)t, U^), (4-30) 

where we have introduced the following functions 

n,{t) = l + t /o (Pe(t')^c - cua) e*'^^*') dt' , v^it) = tX /J e**(*') dt' 

It is not a big surprise that this approximated solution does not share with z{t) and t) 
all their properties. In particular, while for instance [z{t), t)] = for all t, zi{t) and 
Zi(/, t) do not commute. For this reason we consider the equations for n{t) and k{t) as 
in fl4.25p as far as possible, replacing z{t) and Z{f,t) with zi(t) and Zi{f,t) only at the 
last step. 

It is easy to find that the mean values of the first two equations in fl4.25p can be 
written as 

n{t) = '^ = -2XQ{u{zit)Z^a,t))}, 

k{t) = '-^ = 2\p^m{co{zit)z^if,t))}, ^ • ^ 

which in particular implies a well known identity: Pc{t)n{t) + k{t) = for all t, which in 
turns implies that n(t) = Pc{t) n{t). It should be remarked that, because of this relation, 
since M = O implies Pc{t) = -Pc(O) = M, then when M = O the dynamics of the portfolio 
of r is trivial, n(t) = 11(0), even if both n{t) and k{t) may change in time. 

It is now at this stage that we insert zi{t) and Zi{f,t) in the differential equations. If 
uj is the usual number state, and if we call for simplicity 

u{l):=u{zz^ [Z^{J),Zif)]}, 

uj{2):=uj{Z{f)ZHj)[zKz]}, (4.33) 

r{t) = r]i{t) fi2{t) + iv{2) r]2{t) f/i(t) 

then we get 

n-2XQ\ ['r(t')df\ , 

J (4.34) 
k + 2X'^lj^P,{t')r{t') dt'\ . 
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The time dependence of the portfoho can now be written as 

n(t) = n(o) + 5n(t), (4.35) 

with 

6U{t) = n{0 - M) sin2(At) + 

+ (^-2\ $5 1^ r{t') rft'l (M + (O - M) sin2(At)) + 2A 5^ |^ P^t') r{t') dt'^ , (4.36) 

which gives the variation of the portfoho of r in time. We observe that, as it is expected, 
(5n(t) = if A = 0. 

In the last part of this subsection we look for particular solutions of this system under 
special conditions. A more detailed analysis of these results will be discussed in another 
paper, which is now in preparation and where a more general settings will be considered. 
A first remark concerning (14.361) is the following: if O > M, it is more likely for r to have 
a positive 511 (t) if the number of the shares n in his portfolio at time t = is large: if 
at t = the supply of the market is larger than the price of the share then for a trader 
with many shares it is easier to become even richer^. If, on the contrary, O < M, having 
a large number of shares does not automatically produce an increment of the portfolio. 

Coefficients uj{l) and uj(2) can be found explicitly and depend on the initial conditions 
of the market. If, for simplicity's sake, we consider A = {ko}, that is if the reservoir 
consists of just another trader interacting with r, then we get 

^(1) = \f{K)\' (1 + n) k^-''^ {n' - (1 + n') kt''^) 

and 

u{2) = \f{ko)\^ (1 + n') A;{-^^> [n k^+''^ - {1 + n) k^-''^) . 

It is clear that these coefficients coincide ii k = ko and n = n' . 

Let us first fix M = 1, O = 2, A = 1, = c^c = 1, = = 2. Then the plots of 
6Il{t) below, in which n is fixed to be 10, are related to the following different values of 
00(1) and 00(2): {uj{l), uj{2)) = (1, 1), (1, 10), (10, 1). 
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Figure 1: 6U{t) for n = 10 and (w(l),u;(2)) = (1, 1) (left), (w(l),u;(2)) = (1, 10) (middle), 
(u;(l),^(2)) = (10,l) (right) 

The plots do not change much if we fix n = 5 and, surprisingly enough, also the ranges 
of variations of 6Il{t) essentially coincide with those above: n seems to play no crucial 
role here! In Figure 2 we plot 6Il{t) in the same conditions as before, but for n = 5. 

-J V .J V 

Figure 2: 5Il{t) forn = 5 and (cu(l), cj(2)) = (1,1) (left), {uj{l),u{2)) = (1,10) (middle), 
(a;(l),a7(2)) = (10,l) (right) 
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From both these figures we see that, for trader r, the most convenient situation is 
{uj{l),uj(2)) = (1,10): in this case there is only a small range of time in which 6Il(t) 
is negative. For all other times 6Il{t) is positive and n(t) increases its original value. The 
situation is a bit less favorable for other choices of (u;(l), ci;(2)). This is not surprising 
since and uj{2) are related to the initial values of the stock market we are considering 
and, how it is well known, different initial conditions may correspond to quite different 
dynamical behaviors! 
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Now we change the relation between M and O. Therefore we fix M = 2, O = 1, A = 1, 
uja = = 1; = = 2. Again the plots of (5n(t) below are related to the following 
values of: cu(2)) = (1, 1), (1, 10), (10, 1), and we fix n = 10. 



l\ , ... 



' ^ ' r^- — ' — 



A 



I r \ fi I If \ y\: ■\i 4 



Figure 3: 6U{t) for n = 10 (cu(l), cj(2)) = (1,1) (left), (u;(l), cj(2)) = (1,10) (middle), 
(a;(l),cu(2)) = (10,l) (right) 

We see that these plots look very much as those in Figure 1 refiexed with respect to the 
horizontal axis. This means that, for n = 10, the main contribution in 04.361) is the term 
n{0 — M) sin^(At). Of course this is even more evident if n is larger than 10, while for 
small values of n the role of the other contributions in fl4.36p is in general more relevant. 

We have already stressed that, if M = 0, then 511 (t) = for all t > 0. Therefore we 
don't plot (5n(t) in this condition. Instead of this, we finish considering what happens if 
we change the values of uJa and ujc with VLa and VLc- For that we fix, as in the first case, 
M = 1, O = 2, A = 1, while we take uja = = and VLa = = 1- The related plots 
are 




This result is particularly interesting since it shows that, if = uj{2) = 1 and also 
for n small enough, trader r can only improve the value of his portfolio, no matter the 
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Figure 4: 6U{t) for n = 10 cj(2)) = (1,1) (left), (cu(l), cj(2)) = (1,10) (middle), 

(u;(l),^(2)) = (10,l) (right) 

value of t. Remember now that uj{l) = uj{2) is true if the initial conditions for the trader 
T and for the trader of the reservoir, a, coincide. Therefore, this suggests that the relation 
between the parameters uJa,uJc and Qa, is crucial to determine 11 (t) and in particular 
that if we take Qa = > = t is, in a certain sense, in a better condition with 
respect to trader a. It is therefore natural to associate these parameters, for instance, to 
a sort of information reaching the traders, in analogy with the interpretation discussed 
in [1]. We avoid details here since a deeper analysis is needed for a better understanding 
of the role of all the parameters in H. In a forthcoming paper we will focus our interest 
exactly on this point: we will discuss the solution of system fl4.25p under several different 
conditions, and we will also consider different expressions of Pc{t), arising from some 
different economically reasonable hamiltonian or by experimental data. 

V Conclusions and outcome 

In this paper we have carried on the analysis of a stock market in terms of Heisemberg 
dynamics which we began in [1] . In particular, we have generalized the model introduced 
in [1] by introducing a real dynamical behavior for the price of the shares. This is, in our 
opinion, a big achievement with respect to our previous results. 

Section III is just a pedagogical two-traders model which is useful to fix some general 
ideas and giving some definitions. The same model is further generalized in Section 
IV where a non trivial market has been introduced. We have considered two different 
approximations of this model. The first approximation, the so-called stochastic limit 
approach, is useful to get conditions for the staticity of the portfolio of a given trader. 
The second approximation, more useful for the analysis of the general time evolution of 
the portfolio, produces many results and is quite interesting in view of future applications. 

In particular, in a close future we plan to add more kind of shares within the model, 
and to use system (14.22^ with different functions Pc{t), deduced from other hamiltonian 
models or by experimental data. A more long-distance program also includes, for a market 
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with more shares, an analysis of the role of the Heisenberg dynamics in the analysis of 
stock markets. This will be undertaken clearly via a comparison between our results and 
the experimental data. 

We would also like to comment that, as already briefly discussed in [11], the same 
general strategy seems of some utilities in contexts which are apparently very far from 
stock markets and particle physics. Indeed, the mechanism analyzed here is natural 
whenever we are interested in describing exchanges between different active components 
of our (physical, biological, economical,...) system. Indeed, this is exactly the original 
remark which produced second quantization in elementary particle physics, [6]. Just as a 
different example, we may also use the hamiltonian in (14.11) . or some modification of this, 
for a predator-prey system. In this case n represents the number of predator operator 
while k is the number of prey operator. The mechanism in Hj implies that when the 
number of predator increases of one unit the number of preys decreases of < P > units, 
and P can now be interpreted as a sort of ability of the predator to catch its victims. 
We refer to [TT], and to a paper in preparation |T2], for more applications to sociological 
contexts. 
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Appendix A: the definition of 

In this Appendix we will discuss in detail how to define the operators and c^^, and 
some useful formulas related to them. 

To make the situation simpler, we just neglect here the role of the other operators 
appearing in Section III, i.e. the number of share and the supply operators, since they 
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play no role in the definition of, say, . We will just consider two sets of bosonic operators 
c and with [c, c^] = = 1, and the common vacuum vector ifQ-. c(fQ = pipo = 0. In 

a standard fashion we call 

^k,m=^^ic^npT^O, (A.l) 

V k\ ml 

where k,m >0. It is well known, 0, that ipk,m is an eigenstate of /c = c^c and P = p'^p: 
k(pk,m = kipk^m and P (pk,m = f^'^k^m- Siuce we have, if k is large enough, apk,m = 
\/kipk-i,mt c^^k,m = ^/k{k — l)v9fc_2,m, and so on, it is natural to define 

(pk,m, if m = 0, VA; > 0; 

</?fc,,n := <! 0, if m > A;, VA; > 0; (A.2) 

A/A?(A5^n7^^^^TA5^-mrnO '/'fc-m.m, if A; > m > 



Analogously, since cVfc.m = Vk+ l(fk+i,m., (c"^)Vfc,m = V(A; + 1)(A; + 2)ipk+2,m, and 
so on, for all k and m > 0, we put 

c^^ifk ■= I ifm = 0, VA;>0; ^^^^ 

a/(A; + 1)(A; + 2) ■ ■ ■ (A; + m) V5fc+m,m, if m > 0, 

Remark: We could use a different name for the operators and c^^. For instance we 
could call Y and W the operators defined as in (1A.2I) and ( 1A.3I) : however we have decided 
to keep this notation to stress the role of both c and P in the definition of these ladder 
operators. 

These definitions, other than natural, have two nice consequences: (i) they really define 
the operators c and c' since they are now defined on the vectors of an orthonormal 
basis in the Hilbert-Fock space 7i of the system, which is the closure of the linear span 
of the set {'Pk,m, k,m > 0}. In this way we by-pass the problems raised in [1], and we 
can avoid replacing the hamiltonian fl2.1l) with the approximated hamiltonian fl2.5l) : (ii) 
we get an extra bonus which suggests that flA.2p and flA.Sp are good definitions: indeed 
we find that 

and we omit the proof of this claim here. 
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More relevant for us is to deduce some commutation rules which involve the operators 
fc")^, where c" can be c or We claim that 



P,c^ 



t p 

c' c, c 



0, 



-Pc' 



-c^P 



(A.4) 



Again, we omit the proof of these rules here since they can be easily deduced applying 
both sides of each line above to a vector ipk,m of our orthonormal basis. We simply 
remark that, for instance, [A;, c^] = —Pc^ is an extended version of [k, c'] = — / c', while 
[A;,c^'^] = Pc^^ extends [A;,c^'] = lc^\ 



Appendix B: Few results on the stochastic hmit 

In this Appendix we will briefly summarize some of the basic facts and properties con- 
cerning the SLA which are used in Section IV. We refer to [S] and references therein for 
more details. 

Given an open system iS + 7^ we write its hamiltonian H as the sum of two contri- 
butions, the free part Hq and the interaction \Hi. Here A is a coupling constant, Hq 
contains the free evolution of both the system S and the reservoir 7^, while Hi con- 
tains the interaction between S and TZ. Working in the interaction picture, we define 
Hiit) = e'"°*Hie-'^ot g^^^^ (.galled wave operator Ux(t) which is the solution of the 

following differential equation 

dtUxit) = -iXHi{t)Ux{t), (B.l) 

with the initial condition Ux{0) = 1. Using the van- Hove rescaling t ^ see |9] for 
instance, we can rewrite the same equation in a form which is more convenient for our 
perturbative approach, that is 
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with the same initial condition as before. Its integral counterpart is 

which is the starting point for a perturbative expansion, which works in the following way. 

We will limit ourself here to consider the zero temperature situation. Then let ipo 
be the ground vector of the reservoir and C, a generic vector of the system. Now we put 
(y9Q^^ = (y9o ® ^. Wc waut to computc the limit, for A going to 0, of the first non trivial 
order of the mean value of the perturbative expansion of U\{t/X'^) above in (p'^\ that is 
the limit of 

for A — > 0. Under some regularity conditions on the functions which are used to smear 
out the (typically) bosonic fields of the reservoir, this limit is shown to exist for many 
relevant physical models, see [9] and ^0\. We define I{t) = limA-^o -^a (^) • In the same 
sense of the convergence of the (rescaled) wave operator U\{^) (the convergence in the 
sense of correlators), it is possible to check that also the (rescaled) reservoir operators 
converge and define new operators which do not satisfy canonical commutation relations 
but a modified version of these, [10]. Moreover, these limiting operators depend explicitly 
on time and they live in a Hilbert space which is different from the original one. In 
particular, they annihilate a vacuum vector, rjo, which is no longer the original one, ipQ. 

It is not difficult to deduce the form of a time dependent self-adjoint operator Hj^''\t), 
which depends on the system operators and on the limiting operators of the reservoir, 
such that the first non trivial order of the mean value of the expansion of f/^ = 1 — 
i Hj'^''\t')Ut'dt' on the state ?7q^^ = rjQ ® ^ coincides with I{t). The operator Ut defined 
by this integral equation is called again the wave operator. 

The form of the generator follows now from an operation of normal ordering. More in 
details, we start defining the fiux of an observable X = X ® 1,., where is the identity 
of the reservoir and X is an observable of the system, as jt{X) = U^XUf. Then, using 
the equation of motion for Ut and U}, we find that dtjt{X) = iU}[Hj'^''\t), X]Ut. In order 
to compute the mean value of this equation on the state ?7q^'', so to get rid of the reservoir 
operators, it is convenient to compute first the commutation relations between Ut and the 



32 



limiting operators of the reservoir. At this stage the so called time consecutive principle 
is used in a very heavy way to simplify the computation. This principle, which has been 
checked for many classes of physical models, [9], states that, if /5(t) is any of these limiting 
operators of the reservoir, then 

[(3{t),Ut'] = 0, for all t > t' . (B.5) 

Using this principle and recalling that r/o is annihilated by the limiting annihilation oper- 
ators of the reservoir, it is now a simple exercise to compute (dtjAX)) (o and, by means 
of the equation (dtjAX)) (e) = ijALiX))) (o, to identify the form of the generator of the 
physical system, which allows us to obtain equations of motion in general much easier 
than the original ones, since the reservoir disappear. 
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